Application No. 09/786,274 

AMENDMENTS TO THE SPECIFICATION 
In the Specification 

Please substitute the following amended paragraph(s) and/or section(s) (deleted 
matter is shown by strikethrough and added matter is shown by underlining): 
Page 2, line 12-line 27, 

The calculation of the movement vectors of the nodes of the meshwork can be 
carried out according to various methods. First of all there are putting into correspondence or 
"matching" methods consisting of testing a discrete set of possible values of movement vectors 
for each node of the meshwork and of retaining the best vectors according to a given criterion. A 
second method known as a transformed method consists of using the properties of the Fourier 
transform and its extensions so as to convert the movement into a phase jump in the transformed 
space. Finally, there is a third method known as a differential method for determining the 
movement vectors by optimising optimizing a mathematical criterion (for example a quadratic 
error between the image and its aforesaid value with the movement field). This method is most 
frequently used for estimating movement with mod e lisation modelization by finished elements. 
A conventional differential method for optimising optimizing movement vectors is the Gauss- 
Newton method. The present application concerns more particularly the movement estimation 
method family using a model of finished elements and a differential method to determine the 
movement field. 

Page 3, line 15 - page 4, line 17, 

One object of the invention concerns a method for estimating movement in which 
the meshwork is optimis e d optimized during estimation so as to obtain at the end of the method 
of meshwork adapted to the semantic contents of the images. To this effect, the finished 
elements are refined during the movement estimate. 
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* - * Another object of the invention is to improve the effectiveness of the Gauss 

Newton method so as to optimis e optimize the movement vectors of the nodes of the meshwork. 
To this effect, this optimisation optimization is carried out on several resolution levels of the 
images. 

Finally, another aim of the invention is to provide a method for estimating 
movement so as to avoid aforesaid pathological situations. To this effect, the invention provides 
adding during the movement vectors optimisation optimization step constraints so as to avoid 
these situations. 

Also, the invention concerns a method for estimating the movement between two 
numericals image I x and I 2 with luminance Y x and Y 2 and intended to generate for each point of 

— > A 

coordinates x,y of the image h a movement vector d(x,y) = (d xX d y ) so as to form an image 1 2 
from the image l\ with luminance Yi{x,y)-Y x {x - d xX y -d y ) which is an approximation of the 

image h, charact e ris e d characterized in that it comprises the following steps: 

(a) defining an initial model of finished elements comprising a meshing 
whose nodes are points of the image h, a movement vector at each node of said meshing, and an 
interpolation formula calculating the value of the movement vector of each point of the image I2 
from the values of the movement vectors of the nodes of the mesh to which it belongs, 

(b) optimising optimizing the value of the movement vectors of the model 
according to a differential method, 

A 

(c) calculating a variation E between the image h and the image 1 2 for each 
finished element or mesh, 

(d) carrying out a finer meshing on a discrete fraction of the set of finished 
elements determined according to a criterion relating to the variations E and allocating a 
movement vector to each new meshing node, 

(e) repeating the steps (b), (c) and (d) on the model of finished elements 
obtained at the end of the preceding step (d) until a stoppage criterion is satisfied. 
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Page 4, line 25 - line 30, 

Finally, according to a preferred embodiment, constraints are added to the 
movement of the finished elements at the time of optimising optimizing the movement vectors so 
as to avoid the reversal of the finished elements. According to another embodiment, it is also 
possible to introduce constraints so as to avoid the flowing over of the meshing obtained after 
applying the movement vectors beyond the sphere of the image l\ . 

Page 6, line 16 - line 20, 

According to the invention, as the movement is gradually estimated, the value of 
the movement vectors of the meshing nodes known as nodal vectors is optimis e d optimized and 
the meshing is locally densified when this is necessary. Advantageously, this optimisation 
optimization shall be carried out on several resolution levels starting with a low resolution level. 

Page 7, line 13 -line 30, 

According to step (b), the value of the movement vectors of the model are 
optimis e d optimized according to a differential method, such as the Gauss-Newton method or its 
Marquardt extension. This optimisation optimization can be free, that is without constraints 
imposed on the possible values of the nodal vectors, or with constraints. The nodal vectors 
denote the movement vectors of the nodes of the meshing. Optimsation Optimization with 
constraints is directly linked to free optimisation optimization and forms the subject of an 
embodiment given in detail later in the text. 

The free optimisation optimization technique used here makes use of the advantageous 
characteristics of the Gauss-Newton method (rapid convergence of the optimum) and the 
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gradient method with adaptive step (global convergence towards a local optimum) to resolve the 

linear system to be followed. This technique is an iterative correction of the movement vectors 
— » 

d(x, y) making it possible to obtain from the start a rough approximation of the movement. The 
number of iterations k for this optimisation optimization of the movement vectors is either 
specified by the user at the start of the method, or depends on a threshold linked to the maximum 
variation between two consecutive values of nodal vectors for two successive iterations. Below 
are details of the Marquardt extension of the Gauss Newton optimisation optimization method. 

Page 9, line 13 - line 15, 

At the end of this optimisation optimization phase, there are available N nodal 
vectors each relating to one node of the meshing. 

Page 11, line 25 -line 29, 

In practice, the images lj are obtained by the filtering of the image /, using a linear 

low pass filter only allowing Vi of the spectral band of the image in question in the directions x 
and y, that is a pulse response filter h r n having a pass-band BP r = [l / 2 r+1 ,1 / 2 r+1 ) in the space of 

standardis e d standardized frequencies [-1/2, '/?/. The image /f is defined by the following 
equation: 

Page 12, line 10 -line 13, 

As indicated previously, this optimisation optimization on several resolution 
levels is able to improve and accelerate the convergence of the calculations of the movement 
vectors. It is to be noted that the number of resolution levels R selected may differ from the 
number of successive refinings carried out on the meshing. 
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Page 13, line 7 - line 9, 

According to the invention, this constraint is associated with each triangle at the 
time the movement vectors are optimis e d optimized . The step of optimising optimizing the 
movement vectors amounts to a system of the type: 

Page 13, line 15 -line 19, 

So as to resolve the optimisation optimization problems under constraints, the so- 
called increased Lagrangian technique is used. This technique is described in the work entitled 
"Theories and algorithms" by Michel Minoux, Volume 1, pp 257 260 256-259 , published by 
Dunod 1983. This technique combines two optimisation optimization techniques: Lagrangian 
optimisation optimization and the optimisation o ptimization of external penalties. 

Page 14, line 5 - line 6, 

The constraints g e have been previously lin e aris e d linearized by the Taylor 
formula to the order 1 : 

Page 14, line 9 -line 17, 

The optimisation optimization method is then the following: 

- k = 0 is initialis e d initialized 

- X = 0 and r = 0 are placed, X e 9? m and r e 5R m , and m denotes the number of 
triangles of the model 

then the minimum SD k+] {X y r) is determined so that on iteration 
H.5D k+x =V£* -C> 
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where y' = (A,r), C is a matrix of $R 2A x9i 2 
C'y forms a matrix of the lin e aris e d linearized constraints having for coefficients 
the following values: 
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